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ABSTRACT 

We present observational results of the [Fe II] A 1.644 fim emission from the jets 
of L1551 IRS 5. The data sets were obtained through 13 fully sampled slits aimed 
at the base of the jets. These sets are used to construct a three-dimensional cube. 
The field of view was 5"8x4"2. We confirmed that the position of the knot PHK1 
coincides with a stationary, point-like x-ray source within ±(y.'3. The northern and 
southern jets are distinguished from each other at a point 0/6 away from their driving 
sources. We also confirmed that the northern jet consists of well-separated high- and 
low-velocity components (HVC and LVC, respectively). The HVC has a terminal 
velocity of ~400 km s _1 and shows a consistently narrow velocity width of 40 km s _1 . 
The LVC covers the velocity range from Vlsr =0 to —240 km s _1 and has broad 
velocity widths of ~ 150-1 80 km s _1 . These decrease with distance from the driving 
sources. The spatial width of the LVC varies from / /6-0 / /7 at Vlsr ~ —200 km s _1 
to 0/8-0/9 at Vlsr ~ — 30 km s -1 . These characteristics are well understood in terms 
of the two types of outflow mechanisms that are working simultaneously: one is the 
HVC, which is launched in a narrow, inner radial region at 0.04-0.05 AU, and the 
other is the LVC, which is launched in a wider, outer radial region from within 0. 1- 
4.5 AU of the accretion disk. Part of the LVC emission could arise in the gas entrained 
or shocked by the HVC. We also discuss the possibility that part of the HVC gas is 
thermalized at PHK1 to produce the x-ray emission and LVC. 

Subject headings: ISM: Herbig-Haro objects — ISM: individual (HH 154, L1551) — 
ISM: jets and outflows — stars formation — stars: pre-main- sequence 
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1. INTRODUCTION 



L1551 IRS 5 is one of the best objects for the study of outflows and jets emanating from 
low-mass young stars. This is due largely to its proximity (d ~ 140 pc), visibility, and isolation. 
It was the first object in whi ch a molecular molecular outflow , with a CO radial velocity of | V r \ < 
20 km s _1 , was discovered (Snel l. Loren. & P lambeck 1980). T he CO outflow em ission shows an 
apparent linear acceleration ( Uchida et al.lll987 : Moriarty-Schieven & Snellfl988 ). which suggests 
that the outflow should consist of ambient material swept up by an invisible fast jet or wind. A high- 
veloc i ty ( | V m a , Y | ~ 1 00 km s ~ 1 ) at omic wind was later detected inside the CO lobe (ILizano et al 



19881 : iGiovanardi et al.N 19921 |2000j). The existence of a fast wind can be deduced from the many 
optical and H 2 em ission-line features showing large proper motions and radial ve locities inside the 



blue- shifted lobe (IStocke et al.lll988l : iGraham & Heverill990l : 
also the first ob 



Mundt & Fried 



ect to ward which an optical jet was detected (|Cohen. Bieging. & Schwarta ll982; 



Davis et al. 



1995). L1551 IRS 5 is 



1983b . High angular resolution radio observations revealed that L1551 IRS 5 



is a binary, or even triple, system of embedded young stellar objects, referred to as the L1551 



IRS 5 V LA sources, with two jets emanating from them (ILim & T akakuwa 2006: [Rodriguez et al. 
2003aUbl : lltoh et alJEoOol : lHartigan et aDboOol : iFridlund & Liseaulll998l) . 



We show in Figure Q] an image of the jets taken with an [Fe II] narrow-band filter on 2002 
November 27 for reference. It shows the two jets emanating from LI 551 IRS 5 in the [Fe II] 
A 1.644 /im emission. The jets can be distinguished from each other at a position more than 3" 
away from the driving sources in the fi gure. Along th e northern jet, three strong knots (PHK1, 
PHK2, and PHK3) have been labeled (IPyo et al.1 120021) . Several knots in t he northern jet show 
proper motions along the flow, indicating episod i c mass ejection from IRS 5 iFridlund et al. 2005 : 
Fridlund & Liseaull994l : lNeckel & Staudell 19871) . iBachiller. Tafalla. & Cernicharoli 19941) detected 



evidence of episodic mass loss in the molecular outflow along the axis of the blueshifted lobe. The 
molecular outflow has a longer dynamical time scale (a few x 10 4 yr) than the optical jets (10 2 — 10 3 
vr). The episodic mass-loss e vents could be related to the fact that IRS 5 is an FU Ori type star 
iCarr. Harvev. & Lestell987h . 



We have been studying the kinematics and structure of jets at their bases by analyzing the 
near-infrared [Fe II] (A 1.644 /zm) emission line which, due to its small extinction, traces jets 
closer to their driving sources than optical forbidden lines would. This is the third of a series of 
papers on high spectral and spatial resolution [Fe II] spectroscopy of L1551 IRS 5 with SUBARU. 



an 8-meter- diameter ground-based optical and infrared telescope. In previous papers (|Pyo et al. 



2002L |2005|) . we explained that the northern jet has two distinct velocity components: the low- 
velocity component] (LVC) with a velocity of ~ —100 km s _1 and the high- velocity component 



'The LVC in this paper should include the pedestal and red wing component (PWC), identified by 



Pvo et al. 



(2002, 
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(HVC) with a velocity of 300 km s _1 . The LVC is much stronger than the HVC at PHK1, 

the brightest [Fe II] knot located close (1"2 in 2002 December) to the L1551 IRS 5 VLA sources 
(hereafter called as IRS 5 VLA, a singular proper noun). It shows a wide spatial width and a broad 
line width of 130-160 km s _1 near IRS 5 VLA, implying a large opening angle of 50-10 0°. These 
chara cteristics are similar to the molecular Extremely High Velocity (EHV) outflows jBachiller 
19961) . The HVC, on the other hand, consistently shows a narrow line width (~40 km s _1 ), is 
extended farther away from IRS 5 VLA, and is dominated in intensity at PHK2, a bright [Fe II] 
knot that was located ~4'.'7 away from IRS 5 VLA. However, the previous observations did not 
fully cover the [Fe II] emission region, and the morphology and kinematics of the HVC and LVC, 
and their interrelation have remained unclear-especially at PHK1, where the northern and southern 
jets merge together spatially and the LVC is no longer separated from the southern jet because the 
two have a similar velocity. 

The presence of two velocity components in jet emissions is an important clue to under- 
standing the launching mechanism of outflows from young stellar objects. Magnetocentrifugal 
acceleration has been widely accepted as a p romising explanation for the launching mechanism of 



outflows from a young s tar-disk system (see lKonigl & Pudrit z 2000; Shu et al.ll2000l : IPudritz et al 



20071 : IShang et al.ll2007l) . Th ere are two popula r models based on this magn etocentrifugal mech - 
anism: the "X-wind model" (IShang et all 120071) and the "disk-wind model" (IPudritz et al.ll2007l) . 
The X-wind model postulates that an outflow originates from a very narrow radial range of an ac- 
cretion disk near its inner edge, where the highly concentrated stellar magnetic field interacts with 
the disk. The disk-wind model assumes that an outflow is launched from the wider radial range of 
an accretion disk. The terminal velocity of a magneto centrifugally accelerated outflow scales with 
the Keplerian rotation velocity at its launching radius (|Konigl & Pudritzll2000l) . Thus, the velocity 
width of a disk- wind should naturally be wider than that of an X-wind. This difference between 



the tw o models is best recognized in synthesized position-velocity diagrams (PVDs) (|Pyo et al. 



2006). High spatial- and velocity-resolution spectroscopy is essential for the direct comparison of 
observations with these model calculations. 

Here we present the results of [Fe II] line slit-scan observations covering the base of the LI 55 1 
IRS 5 jets with 13 fully sampled transverse slit locations. We constructed a three-dimensional data 
cube to analyze the spatial-velocity structure of the [Fe II] jets. 



2005), of the northern jet based on the results of this paper. 
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2. OBSERVATIONS AND DATA REDUCTION 



Observations of the L1551 IRS 5 jets were conducted on 2004 December 1 (UT) using the 
Infrared Camera and Spectrograph (IRCS. iKobayashi et alj|2000c iTokunaga et al.ll 19981) mounted 
on the Cassegrain focus of the Subaru Telescope atop Mauna Kea, Hawai'i. The //-band spectra 
were taken using the IRCS echelle spectrograph, which is equipped with a Raytheon 1024x 1024 
InSb array with an Aladdin II multiplexer. The near-infrared slit viewer of IRCS allowed us to 
place the slit accurately at selected positions on the jets. This was essential to match spectra with 
spatial coordinates. We used a slit of / /6x5'. / 9 (width x length). The slit width provided a velocity 
resolution of 60 km s _1 (A/AA =5,000). The pixel scale along the slit length was / /060 pixel -1 . 
The spectra were taken around PHK1 at 13 slit positions, at a position angle (PA) of 346°, which 
was perpendicular to the northern jet. Adjacent slit positions were separated by 0"3, i.e. half of 
the slit width. The field of view covered by the 13 slits was 5! 9 x 4.'2. We made three sets of slit 
scans. Because of the variable sky conditions, only one set showed good data quality. That set was 
used for this paper. The total on-source exposure time was 180 s for each slit position. The seeing 
(FWHM) was ~0 / /4 at the //-band. The area covered by the current slit scan observations is shown 
in Figured] (lower panel) by a thick rectangle, together with the coordinate system (X, Y) that we 
use in this paper. These are all superimposed on the contour versions of the upper panel image. 
The positions of IRS 5 VLA, PHK1, PHK2 and PHK3 are also marked. 

Sky spectra were obtained before and after each set of slit scans. The atmospheric OH airglow 
emission was satisfactorily removed from the object spectra when the properly scaled sky spectra 
above the background level were subtracted from them. In order to find a proper scaling factor for 
each object's spectrum when removing OH features, we used a strong OH line at A 1.6502 /im to 
derive a scaling factor so that the standard deviation of the residual (Object — scaled Sky) at the 
line was minimized, while keeping the average of the residual at zero. We observed the standard 
star HD 35656 (A0V, T c g = 9480 K, V — 6.403 mag) to use its spectrum for removing telluric 
lines and to calibrat e flux. The dat a were reduced mainly using the ECHELLE and LONGSLIT 
packages in IRAF (|Pyo et al.ll2002|) . After making a three-dimensional data cube with IRAF, we 
used IDL to analyze it. 



3. RESULTS 

3.1. Continuum- Subtracted [Fe II] Image 

Figure |2h shows the continuum- subtracted [Fe II] line image integrated over the entire emis- 
sion velocity of —420 to +20 km s _1 with respect to the local standard of rest. Note that this is a 
pure [Fe II] line image in which the continuum emission is subtracted accurately down to the noise 
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level. Slit-scan imaging allowed such accurate continuum subtraction because the continuum and 
[Fe II] line were observed at exactly the same time. 

It is difficult, on the other hand, to construct a pure line image from two narrow-band filter 
images "on" and "off" a given line emission if they are not taken simultaneously. This is because 
the continuous variation of sky conditions makes it difficult to match the continuum levels between 
the two images. 

In Figure^, the [Fe II] emission peak is called PHK1 (Paper I) and is plotted at X =0" and 
Y = 0", where X is along the slit length (PA = 346°) and Y is perpendicular to it (PA = 256°). A 
new knot is seen in Figure [2b at Y = 0?6 near the tip of the elongation. The faint southern jet is 
seen at X ~ — 0/6 and Y > 1". The northern and southern jets are merged and indistinguishable 
at PHK1. The PAs of the northern and southern jets are 261?6 ± 1?2 and 235?2 ± 1?4, respectively. 

Figure [2b shows the continuum emission averaged over the velocity ranges from —1840 to 
—560 km s _1 and from +180 to +1500 km s _1 with respect to the [Fe II] line rest wavelength 
in the LSR. The continuum peak is located at (X, Y) = (0 / /06, —0/6): the offset between the 
continuum peak and PHK1 is 0'.'6 (Act = -0'.'57, A5 = -0'.'20). 



3.2. Velocity Structure 

Figure [3] shows the velocity channel maps of the [Fe II] A 1.644 fim emission for the velocity 
range from Vlsr = —414.6 to +28.0 km s _1 . Each panel gives an integrated intensity map over 
a width of 29.5 km s _1 . The spatial width across the X-axis at the emission peak at each velocity 
channel is listed in Table CD 

The HVC is seen in the panels from[3b to [3b (-381.5 km s" 1 < Vlsr < -267.1 km s" 1 ), 
where the emission is elongated along the F-axis with a spatial width of 0/56-0/68 (FWHM 
not deconvolved with the seeing size) centered around X = ff!15 (see Table LB. The emission 
peak moves from larger to smaller positions of Y, with the velocity in each panel reflecting the 
velocity gradient in the HVC. Although the channel maps in panels [3]/ and[3b (—267.1 km s _1 < 
Vlsr < —208.1 km s _1 ) show an elongated geometry similar to that of the HVC, both the HVC 
and LVC of the northern jet contribute to the emission in these channels (see Figure HI). The 
emission in panel [3pi (—208.1 km s _1 < Vlsr < —178.5 km s _1 ) comes from the LVC without 
contamination from the southern jet. The spatial width at this velocity is 0/66, which is similar to 
that of the HVC. In the panels [3H3 (-178.5 km s -1 < Vlsr < -60.5 km s" 1 ), the southern jet 
overlaps with the LVC and the spatial width becomes as large as 1". Channel maps[3}?i and[3h 
(—60.5 km s _1 < Vlsr < —1-5 km s" 1 ) are free from the southern jet emission and show the less- 
blueshifted, lower- velocity emission from the northern jet. The spatial widths at these velocities 
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are (X'86 and (X'83, respectively, and are significantly larger than those of the LVC in panels |3^ and 
which show more-blueshifted, higher-velocity emission from the northern jet. 

Figure© shows a position- velocity diagram (PVD) where the [Fe II] flux integrated over the 
X-axis (perpendicular to the northern jet) at each position of Y is plotted against the radial velocity 
(Vlsr). This shows the well-separated high- and low-velocity components discussed above: the 
HVC at Vlsr ~ -300 km s" 1 and the low-LVC covering Vlsr from -240 km s _1 to km s~\ 

Figure HJ: shows an X-velocity diagram (XVD) where the [Fe II] flux integrated over the Y- 
axis (parallel to the northern jet) at each position of X is plotted against the radial velocity. An 
emission feature is seen at ~ —120 km s _1 corresponding to the southern jet (marked as "A") in 
this figure. The HVC is also seen at ~ —300 km s _1 around X ~ 0/2 (marked as "B"). 

Figure |4b can be deconstructed along the F-axis to produce Figure [51 a series of XVDs at 
various positions of Y. This is one of the most remarkable diagrams created from the slit scan 
observational data. First, the southern jet can be traced down to Y=— 0/60, i.e. as close as 0/6 
from the driving source IRS 5 VLA. This is because the radial velocity of the southern jet is almost 
constant at Vlsr ~ — 120±10 km s _1 with a relatively narrow FWHM width of 80.5±15.3 km s _1 
(FWHMfl in contrast to the large velocity width of the LVC (150-180 km s" 1 FWHM). 

Even at Y <0?!6 where the LVC mingles with the southern jet, the southern jet can be sepa- 
rated from the LVC. Note that the southern jet can be traced only at Y >1" (>2" from the driving 
source) in the velocity-integrated spatial map in Figure [2fr (or Figure Sh). Figure |5] shows that the 
southern jet gradually drifts away from X = —0/1 to — C/8 as Y increases from —0/60 to l'/62. 
The flux density of the southern jet reaches its maximum at Y=0", i.e. at PHK1, where it is 2.5 
times larger than that at F=l'/62. This means that the southern jet significantly contributes to the 
[Fell] flux ofPHKl. 

Second, the lower (less-blueshifted side) velocity cutoff of the LVC tends to decrease with 
Y. This is seen in Figure [5] as the lowest contours intersect the X= line at Vlsr ~ km s -1 for 
Y < 0"; this occurs at ~ —50 km s _1 for Y > l'/2. The higher (more-blueshifted side) velocity 
cutoff of the LVC is rather constant in this range of Y. An interesting point is that the decrease in 
the lower velocity cutoff of the LVC is similar to that seen in the HVC: its peak velocity decreases 
from Vlsr = —263 to —312 km s _1 with increasing Y from 0" to 1". However, we conclude that 
the two similar tendencies are unrelated because the velocity gradient of HVC seen for 0" < Y <1" 
is time variable (see 513.31 ), while the lower cutoff velocity gradient of LVC does not vary with time, 
as seen in the PVD taken four years ago. 



2 The FWHM width is 53 km s 1 if deconvolved with the velocity resolution of 60 km s 1 
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3.3. Comparison with Previous Data 

In this section we summarize some of the important points for the discussion below by com- 
paring the current [Fe II] data with previous spectroscopic data taken in 2000 December (Paper I) 
and 2002 November (Paper II). 

1. The pure [Fe II] line image (Figure [2h) bears a close resemblance to Figure 2a (Paper I) that 
was constructed from two narrow-band filter images "on" and "off" of the [Fe II] line. This 
means that the continuum was satisfactorily subtracted in the previous work. Because of 
the slightly improved seeing, or possible time variation between 2000 and 2004, Figure \2)a 
shows the elongation of PHK1 along the northern jet more clearly than Figure 2a in Paper I. 

2. The offset between the //-band continuum peak and PHK1 remains the same as before. It 
was 0/56 in 2000 December (Paper l|] and O^O in 2004 December (present work). The 
difference in the offsets is within the measurement error of ~0/l and is much smaller than 



that expected from the proper motions of 0. 2-0. 44 yr 1 for the other knots (|Fridlund et al 



20051 : Fridlund & Liseaulll994t LNeckel & StaudeHl987|) . Because the continuum peak posi 



tion seems stationary, or at least its proper motion is much smaller than those of the other 
knots, PHK1 may also be stationary with respect to the IRS 5 VLA. Previous measurements 
(noted in Paper I) revealed that, within an uncertainty of 0/25, the //-band continuum peak 
coincides with the ^-band peak positio n, which is located at 0/57 ± 0/35 (Aa = —0/55, 
A5 = -0/13) away from IRS 5 VLA dcampbell etalJll988h . This places IRS 5 VLA at 
Y = -L'2±0M. 

3. The Y-velocity diagram (YVD) in Figure Hp is compared with Fig. 2b in Paper I, which 
is a YVD taken in 2000 December, but is not integrated over X at each position of Y. 
Although both YVDs show the two velocity components, there is a notable difference in the 
radial velocities of the HVC. The radial velocity of HVC does not vary largely with Y for 
0" < Y <\" in the 2000 December data (Figure 2b in Paper I), while it decreases from 
-263 to -312 km s" 1 with Y for 0" <Y <\" in the 2004 December data (Figure©). This 
difference is a result of time variation because the radial velocity gradient of HVC seen in 
Paper I is not reproducible from any cross-section of the current data cube (see Figure [3]). 

4. Figure (XVD) is compared with Fig. la of Paper II, which is an XVD taken in 2002 
November but is not integrated over Y. The two diagrams share common features, except 



3 Note that in lPvo et al] (12002b the separation between the continuum peak and PHK1 was reported to be Of! 7, which 
was measured using the inco rrect pixel scale of 0'.'075 pixel -1 . The actual pixel scale was 80% of this value and the 



true separation was / /56 (see lPyo et al 1 l2003h . 
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that the southern jet at Vlsr ~ — 120 km s _1 for X <0" is more conspicuous in Figure 
because it is integrated over Y where the southern jet gives significant contribution. The 
lowest contours in these figures show a very large spatial width corresponding to the spatially 
wide subcomponent of LVC (LVCwide) reported in Paper II. 

5. Figure |5] suggests that the broad, "flat topped" [Fe II] line profiles in LVC that were reported 
in Paper I at Y — — OCA to —0/8 are most likely formed as a result of the overlapping 
emission from the LVC and the southern jet. 

6. Figure [5] shows that the less-blueshifted side cutoff velocity of the LVC decreases (becomes 
more blueshifted) with Y. As a result, the velocity width of the LVC decreases with distance 
from IRS 5 VLA. It is also clear from Figure[5]that this tendency is not significantly affected 
by the southern jet emission. Thus, the velocity narrowing of the LVC along the outflow, 
which was pointed out in Paper I (Figure 2b there), may be related only to the northern 
jet. We interpreted this tendency as due to the collimation of the LVC in Paper I, and that 
interpretation remains valid. 



4. DISCUSSION 

4.1. The Two Velocity Components: HVC and LVC 

The present slit scan observations provided unprecedented spatio-velocity information for the 
two jets emanating from IRS 5 VLA at the spatial and velocity resolutions of OCA and 60 km s _1 , 
respectively. This has enabled us to separate the northern and southern jets even in the close 
vicinity (<1") of IRS 5 VLA (Figure |5]), although the two jets are discernible only at >2" away 
from IRS 5 VLA if no velocity information is available. 

There is no clear evidence to show which of the two IRS 5 VLA sources is responsible for 
the no rthern or southern jets. The northern and southern radio jets detected within ~1" of IRS 5 



VLA (IRodriguez et al.ll2003bl) have PAs of 247°and 23 5 Respectively, which are similar to those 
of the northern (262°) and southern (235°) [Fe II] jets. Although the PA of the northern radio jet 
is somewhat smaller than that of the [Fe II] jet, it bends toward the west at 0/7 west-southwest 
of the northern driving source and its PA changes to ~260°, which agrees well with that of the 
northern [Fe II] jet. If we assume that the northern and southern radio jets are smoothly connected 
to their respective [Fe II] counterparts at ~1" west of IRS 5 VLA, it is natural that the northern 
and southern radio jets are the inner parts of northern and southern [Fe II] jets, respectively. 



While the southern jet shows a single velocity component with a relatively narrow line width 
and little radial velocity variation, the spatial and velocity structure of the northern jet is more 



-9- 



complicated. Below, we summarize the characteristics of the northern jet revealed by our current 
and previous observations of the [Fe II] A 1.644 /mi line. 



1 . The northern jet covers a wide velocity range of —380 to km s 1 . 

2. The position-velocity diagrams along and across the northern jet show two distinct compo- 
nents with different radial velocities: one is the HVC at a velocity of ~ —300 km s _1 and 
the other is the LVC covering the velocity range of —240 to km s _1 . 

3. The HVC shows a narrow line width with a deconvolved FWHM width of ~40 km s _1 (see 
Paper I). 

4. The LVC shows a broad line width of 130-160 km s" 1 (FWHM, see Paper I). This width 
is intrinsic to the LVC of the northern jet and is not caused by the contamination from the 
southern jet (see item 6 of §3.31) . 

5. The LVC shows a spatial width of ~0.'65 at its more-blueshifted velocity side (Vlsr = —200 
km s _1 ). This width is similar to that of the HVC. The LVC has a larger spatial width (~0 / /85) 
at its less-blueshifted velocity side (Vlsr = —30 km s _1 ). 



An important issue to consider is whether the HVC and LVC of the northern jet represent 
two types of jets from a single driving source or merely an inner and an outer part of a single 
jet. We discussed this in Papers I and II, with preference to the presence of two outflows from a 
single driving source. This is because the HVC and LVC are identified as two distinct peaks in 
the PVDs and show different kinematical characteristics as typically seen in their velocity widths. 
The HVC and LVC are, on the other hand, not easily discernible in the channel maps of Figure [3j 
which might suggest that the HVC and LVC are smoothly connected. However, we should note 
that it is difficult to theoretically simulate PVDs with two distinct peaks under the assumption of 
monotonically increasing or decreasing variation of physical conditions from the jet axis to its 
outer annuli. The two distinct peaks in the PVDs may thus imply that there are two different radial 
regions flowing at different velocities in the northern jet, if axial symmetry holds, and that the two 
annuli emit equally strong [Fe II] emission, even if they are part of a single outflow. 

Another point that has to be stressed is that the HVC always shows a narrow line width, while 
the LVC has a larger line width for the four [Fe II] jets so far observed at sufficiently high spatial 
and spectroscopic resolutions. Table |2] summarizes the deconvolved velocity widths (FWHM) of 
the HVCs and LVCs for those [Fe II] jets. The HVCs of the jets from DG Tau, RW Aur, and HL Tau 
are conspicuous and are easily identified in their PVDs. They all show a width of ~50 km s _1 with 
their blueshifted radial velocities of ~200 km s _1 . The L VC of the DG T au jet is extended and its 



broad velocity width is clearly seen in Figures 1 and 2 of [Pyo et al.1 (|2003|) . The LVCs of RW Aur 
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and HL Tau are compact and are detected as less-blueshifted wing or shoulder emission at or near 
the positions of their driving sources. These two LVCs do not show distinct peaks in their PVDs, 
although the L1551 IRS 5 northern and DG Tau jets do. The LVC emission from all these sources 
fills the radial velocity range between the HVC and ~0 km s -1 , suggesting a velocity width of 
100 km s _1 or larger. The southern jet of L1551 IRS 5 is, in this sense, categorized as an HVC. 
This is based on its typical deconvolved velocity width of 53 km s _1 and extended geometry, 
and is as listed in Table |2] The current paradigm of magnetocentrifugally accelerated outflow 
models predicts that the terminal velocity of a jet is m ainly determined by the Alfven radius , 



within which the magnetic field lines work as lever arms (Konigl & PudritzH2000t IShu et al.ll200C : 



Shang. Shu. & Glassgoldlll998l : iKudoh. Matsumoto. & Shibatalll998l) . Thus, the narrow velocity 



width of a jet indicates that the line-emitting gas is launched within a narrow radial range of an 
accretion disk, while the broad velocity width of a jet (after it becomes collima ted) implies a wider 
range of launching radius. On the basis of this assumption, Pyo et ali (|2006|) suggested that the 



HVCs of RW Aur and HL Tau should be launched from a narrow range of radius or that only such 
outflowing gas originating from a narrow radial range should be sufficiently excited to contribute 
to the line emission. The same interpretation can be applied for the HVCs from L1551 IRS 5 and 
DG Tau, judging by the similarities of their HVCs to those of RW Aur and HL Tau. 

On the other hand, the LVCs show a velocity width 2-4 times larger than the HVCs (Table|2]), 
meaning that the LVC-emitting gas is launched from the wider outer radial range of an accretion 
disk if the streamlines are collimated; i.e. the LVC is extended more than a few tens of AU from 
their driving sources. Because this length scale corresponds to ~Qf!2 at the assumed distance of 
140 pc, the spatially extended LVCs observed for the northern jet of L1551 IRS 5 and the DG Tau 
jet are consistent with the outflow launched from a larger radial range of an accretion disk. There 
is also a possibility that the spatially compact LVCs toward RW Aur and HL Tau may have arisen 
in a region where the streamlines are not yet collimated. 

In summary, the presence of spatially extended HVCs and LVCs with very different kine- 
matical characteristics suggests two types of outflows from a single driving source: the HVC is 
launched from an inner, narrow region of an accretion disk and the LVC from an outer, wider re- 
gion. Figure |6] shows a schematic diagram showing the launching radii of the HVC and LVC for 
the northern jet of LI 551 IRS 5. We simply assumed a Keplerian-rotating accretion disk around 
a 1 M star with the Alfven radius three times that of the launching point radius along each mag- 
netic line of force. The HVC with a narrow velocity width is launched in a very small radial range 
from 0.04 to 0.05 AU from the star, while the LVC with a broader velocity width is launched in 
a radial range from 0.1 to 4.5 AU, 400 times larger radial range than that of the HVC. The re- 
gions responsible for the HVC and LVC are spatially separated on the accretion disk. If both the 
HVC and LVC arise from a single outflow launched continuously over a radial range from 0.04 
to 4.5 AU, we must seek physical and chemical conditions on the disk surface to heat up the gas 
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along the two separated groups of streamlines. Instead, if we assume two mechanisms working 
simultaneously-for example, if the HVC is driven by an X-wind type mechanism and the LVC is 
driven by a disk- wind type mechanism-it may be possible to have the two emitting gases along the 
two independent streamline groups as shown in Figure [6] 

One of the interesting results of our current observations is shown in Figure [3l where the 
LVC itself shows a variation in spatial width that varies with radial velocity, as quantitatively mea- 
sured and listed in Table [TJ The spatial width of the LVC is ~0 / /65 in its very-blueshifted velocity 
channels (Figures (3]? and[3]V), while it is ~0 / /85 at its less-blueshifted velocity channels, shown in 
Figures [3}?2 and(3h- The difference in the spatial width of the LVC with velocity provides obser- 
vational evidence that slower- velocity gas flows along outer and more widely-opened streamlines. 
This may support the idea that the LVC (but not the HVC) of the LI 551 IRS 5 northern jet is an 
outflow launched over a wide radial range of an accretion disk. Alternatively, part of the LVC- 
emitting gas flowing in the vicinity of the HVC may arise in the shocked gas entrained by the 
HVC. 

The similar spatial width of the HVC compared to that of the very-blueshifted part (Vlsr ~ 
—200 km s _1 ) of the LVC may imply that the gas responsible for this part of the LVC flows 
immediately outside the HVC with a velocity difference of ~100 km s _1 . The two flows may 
interact to produce the fainter "pedestal" emission connecting the HVC and LVC in the PVDs. We 
discussed in Paper I a possibility that the pedestal emission and the "red wing emission" at Y > 2'!5 
are produced as a result of turbulent entrainment of the LVC gas. The pedestal emission has a 
velocity (Vlsr ~ —200 km s _1 ) similar to the secondary peak velocity (Vlsr ~ — 180 km s _1 ) 



of the Ha emission taken toward the base of the northern jet dStocke et al.l 11988) and the red 
wing emission shows a tendency of apparent acceleration with distance from IRS 5 VLA. These 
characteristics favor the idea of pedestal and red wing emissions arising in the gas entrained by the 
HVC flow. On the other hand, the main part of the LVC does not exhibit such apparent acceleration 
along the flow, although it does show velocity narrowing (Figure [5]), and is not affected by the time 
variation that was observed in the HVC at 0" < Y <\". This suggests that it is difficult to consider 
the entire LVC emission as arising from the entrained gas. 

We also pointed out in Paper II that the spatially wide subcomponent of the LVC (LVCwide) 
cannot be considered as arising from a bow shock because there is no systematic velocity difference 
or gradient between the LVC and LVCwide- On the other hand, the H 2 v=l-0 S(l) emission is 
less blueshifted (Vls r ~ 60 km s" 1 ) than [Fe II] and has a steadily increasing velocity component 



along the flow (IDavis et al. 2001), suggesting that it originates from the gas entrained by the LVC 
or LVCwide- Another possibility that the significant part of the LVC flow is a shock-decelerated 
HVC gas will be discussed below in relation to the x-ray emission. 

Finally, we should note that there is a third object recently discovered at 0/09 southeast of 
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the northern binary component (|Lim & Takakuwal I2006Q . The possibility that this third object 
contributes to either of the HVC or LVC is small because there is no hard evidence of an associated 
radio jet; the HVC and LVC of the LI 551 IRS 5 northern jet are driven by a single source. The 
HVC and LVC show exactly the same PAs as seen in Figure [3] and both are spatially confined 
around the F-axis (X =0") even at Y ~4" (see Figure 3a of Paper II). This level of alignment 
between the HVC and LVC is rather unlikely if different driving sources are responsible for the 
HVC and LVC. The peak of the HVC emission is systematically shifted to the positive X direction 
by (X'l— 072 with respect to that of the LVC at each position of Y. The reason for this is unknown, 
but such a systematical shift may also be unnatural if the HVC and LVC originate from different 
driving sources. 



4.2. The X-ray Source 



The soft x-ray emission associated with LI 55 1 IRS 5 was first detected by iFavata et al.1 (120021) 
with XMM-Newton, although the spatial resol ution was not sufficient to pinpoint th e actual location 
of the x-ray source on the jets. Subsequently. Bally. Feigelson. & Reipurthl (120031) located an x-ray 
source at O.'5-l" west of IRS 5 VLA from Chandra observations. This raised t he possibility that 
the x-ray emission comes from PHK1, and lBally. Feigelson. & Reipurthl (|2003[) proposed several 
models on this basis. Because the actual location of the x-ray source is still uncertain, we re- 
examined the location of the x-ray emission peak. 

We confirm in this paper that the //-band continuum peak, which agrees with the A^-band peak 

within an uncertainty of 0/25, is located 076 east-northeast of PHK1 and the di stance did not change 

betw een the two observing epochs separated by four years. On the other hand. lBally. Feigelson. & Reipurth 
(I2003|) reported that the published infrared position for IRS 5 is about 0. 5 or more east of the x-ray 
source. These two measurements suggest that the peak position of the x-ray emission coincides 
with PHK1 within an uncertainty of ~073. Note that this result does not rely on the radio position 
of IRS 5 VLA, which has a larger positional uncertainty as a result of a larger registration error 
between the optical and radio reference frames. 



Favata et al.1 (|2006Q reported (based on their new Chandra observations) that a stronger, point- 



like component of the x-ray source did not show any detectable proper motion, while a weaker 
component expanded in size by ~300 AU over four years, which corresponds to the projected 
velocity of ~330 kms -1 . It is interesting to compare these facts with another set of findings 
showing that the proper motion of PHK1 was not detected between the two epochs of observations 
in 2000 and 2004, while the HVC of the northern jet has a radial velocity of ~300 kms -1 , which 
is equal to the projected velocity if we assume a jet inclination angle of 45°. Such similarities 
strengthen the possibility that the point-like x-ray emission originates in PHK1 and further implies 
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that the HVC, but not the LVC, is responsible (at least) for the extended x-ray component. The 
extended and rapidly expanding x-ray component might be related to the newly seen knot at F=0/6 
in Figure [21 

The mechanism of the x-ray emission near PHK1 is still unclear. The radial velocity of the 
HVC (260-350 km s -1 ) corresponds to the actual velocity of 370-500 km s -1 with the above as- 
sumed inclination. Such an outflow has sufficiently high energy to pr oduce a plasma of 10 6 -10 7 K 



responsible for the x-ra y emission when thermalized in a fast shock (Bal ly. Feigelson. & Reipurth 



20031 : iRaga et all 120021) . Because the highly collimated HVC is present at PHK1 and in the re- 
gion farther away from the driving source, it may be difficult to assume that the HVC is fully 
thermalized at PHK1. However, the spatial coincidence of the HVC with the LVC (with its large 
velocity dispersion of 240 km s _1 around PHK1) might imply that part of the HVC is thermalized 
to produce the x-ray emission and LVC. In such a case, a significant part of the LVC might be 
shock-decelerated HVC gas. 



5. SUMMARY 

The main results of this paper are summarized as follows: 

1. The LI 551 IRS 5 VLA sources and the //-band continuum peak are located l'/2 and O^, 
respectively, east-northeast of the strongest [Fe II] knot PHK1 along the northern jet. The 
relative position of PHK1 with respect to the //-band continuum peak did not show any 
variation between the two observing epochs in 2000 and 2004, suggesting that PHK1 is 
stationary. 

2. The southern jet significantly contributes to the [Fe II] flux at PHK1. It has a constant 
velocity of Vlsr=— 120 km s _1 with a comparatively narrow velocity width of 53 km s _1 
(deconvolved FWHM), and is traced as close as 0/6 from the IRS 5 VLA sources. 

3. The high- and low- velocity components (HVC and LVC) of the northern jet are clearly sep- 
arated in position- velocity diagrams. The LVC shows a wide velocity coverage from —240 
to km s _1 with a deconvolved FWHM width of 150-180 km s _1 . It has a spatial width of 
/ . / 6-0 / /7, similar to that of the HVC, at its more-blueshifted velocity (~ —200 km s _1 ), while 
it shows a wider spatial width of 0/8-0/9 at its less-blueshifted velocity (~ —30 km s _1 ) 

4. The high velocity and narrow line width of the HVC implies that its launching point is 
located over a narrow radial range in the inner part of an accretion disk, while the wide line 
width of the LVC means that it is launched from a wider radial range if it is also launched 
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from the accretion disk. The variation of the spatial width with velocity observed for the 
LVC suggests that slower gas flows in an outer annulus with a wider opening angle. 

5. The velocity of the HVC shows a time variation between the two observing epochs in 2000 
and 2004. 

6. The position of the stationary, point-like x-ray source with respect to the near-infrared po- 
sition of L1551 IRS 5 indicates that it coincides with PHK1 within an uncertainty of 0"3. 
The projected expanding velocity of the weaker, extended component of the x-ray emission 
agrees well with the HVC velocity. 

We thank our referee, Dr. Rafael Bachiller, for providing us with valuable comments that 
have greatly improved this paper. We are grateful to the entire staff of the Subaru Telescope for 
their dedicated support to the telescope and observatory operations. This research has been made 
using NASA's Astrophysics Data System and the SIMBAD database operated at CDS, Strasbourg, 
France. 
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Table 1 . Spatial width of the northern jet as a function of radial velocity 



Channel 




Peak Flux 


PeakX 


PeakF 


FWHM Width a 


in Figure [3] 


(kms- 1 ) 


(x 10~ 18 Wm- 2 A" 1 ) 


(arcsec) 


(arcsec) 


(arcsec) 


c 


-355.6 to -326.1 


0.13 


0.18 


1.20 


0.58 


d 


-326.1 to -296.6 


0.23 


0.12 


0.90 


0.56 


e 


-296.6 to -267.1 


0.17 


0.06 


0.60 


0.68 


f 


-267.1 to -237.6 


0.12 


0.06 


0.12 


0.77 


9 


-237.6 to -208.1 


0.12 


0.18 


0.84 


0.64 


h 


-208.1 to -178.5 


0.17 


0.12 


0.84 


0.66 


i 


-178.5 to -149.0 


0.22 


0.00 


-0.06 


0.96 


J 


-149.0 to -119.5 


0.38 


-0.12 


0.12 


1.05 


k 


-119.5 to -90.0 


0.44 


-0.18 


0.00 


0.97 


I 


-90.0 to -60.5 


0.38 


0.00 


-0.06 


0.85 


rn 


-60.5 to -31.0 


0.30 


-0.06 


-0.24 


0.86 


n 


-31.0 to -1.5 


0.13 


0.06 


-0.42 


0.83 



Not deconvolved with the seeing size of / /4 
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Table 2. Deconvolved velocity widths (FWHM) of the HVCs and LVCs for [Fe II] Jets 



HVC LVC 



Jet 


Width 


Remark 


Width 


Remark 


Reference 




(km s- 1 ) 




(kms- 1 ) 






L1551 IRS 5 Northern Jet 


40 


Extended 


150-180 


Extended 


Present work 


DGTau 


50 


Extended 


~100 


Extended 


Pvo et al. (2003} 


RW Aur 


50 


Extended 


100 


Compact 


Pvo et al. (2006) 


HLTau 


40 


Extended 


>ioo 


Compact 


Pvo et al. (2006) 


LI 551 IRS 5 Southern Jet 


53 


Extended 


N/A 


Not Detected 


Present work 



Fe ll] + cont 



2 




6RA (") 



Fig. 1. — (upper) An [Fe II] narrow-band filter (A c =1.644 fim, AA =0.026 /mi) image of the 
L1551 IRS 5 jets taken on 2002 November 27. Note that the emission includes both continuum 
and [Fe II] line, (lower) The area covered by the present slit-scan observations is shown with a 
thick solid-line rectangle superimposed on the contour version of the upper image. The continuum 
peak, PHK1, PHK2 and PHK3 (Knot D) are indicated. The red dotted-lines trace the brightest 
parts of the no rthern and southern je ts. The thick plus sign marks the mid point of the L1551 IRS 5 
VLA sources ([Campbell et al.lll9 88). The North is up and the East is left. 
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Fig. 2. — (a) Continuum-subtracted [Fe II] line image integrated from —420 to +20 km s -1 in 
Vlsr- The upper panel shows the directions of the northern and southern jets on a. The dash-dotted 
horizontal and vertical lines define the peak of PHK1 at (0", 0"). The thick plus sign marks the 
midpoint of the L1551 IRS 5 VLA sources Jcampbell et alJll988h . The contours are drawn with 
equal intervals in a logarithmic scale. The contour levels are 0.12976, 0.30828, 0.51061, 0.73992, 
0.99980, 1.29435, 1.62817, 2.00650, and 2.43529 x 10~ 18 W m~ 2 . The filled circle at the lower 
right corner shows the seeing size of 0. 4. (b) Continuum emission averaged over the two line-free 
channels of —1840 to —560 km s _1 and +180 to +1500 km s _1 . The dash-dotted horizontal and 
vertical lines define the peak of the continuum at (0/06, —0/6). Other signs and contour levels are 
the same as in a. 
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Fig. 3. — Channel maps of the [Fe II] A 1.644 /zm emission. The velocity interval of each 
channel map is 29.5 km s _1 . The contours are drawn with equal intervals in a logarithmic scale. 
The contour levels are 0.01420, 0.05368, 0.09472, 0.13737, 0.18170, 0.22779, 0.27569, 0.32549, 
0.37725, and 0.43106 x 10" 18 W m~ 2 . The filled circle at the lower right corner of panel a shows 
the spatial resolution (seeing size) of 0"4. 
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Fig. 4. — The data cube is shown from three points of view, (a) Same as Figure [2fr. (b) The Y- 
velocity diagram integrated over the X-axis at each Y. The contour levels are 0.23560, 0.40585, 
0.59604, 0.80863, 1.04626, 1.31187, 1.60877, and 1.94064 x 10~ 18 W m~ 2 A -1 , (c) The X- 
velocity diagram integrated over the F-axis at each X. The contour levels are 0.23620, 0.55682, 
0.93978, 1.39721, 1.94360, 2.59623, 3.37578, and 4.30693 x 10" 18 W m" 2 A -1 . Label A shows 
the southern jet emission (Vlsr = — 120 km s _1 ) and Label B marks the high velocity component 
(Vlsr = —300 km s _1 ) of the northern jet. The contours are drawn with equal intervals in a 
logarithmic scale. The filled ellipse at the lower right corner of each panel indicates the resolutions 
of 0. 4 in space and 60 km s _1 in velocity. 
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Fig. 5.— X-velocity diagrams at Y = (a) -l'.'l, (b) -0"6, (c) 0"0, (d) 0"6, (e) l"2, and (f) l"62. 
The contour levels are 0.00627, 0.01667, 0.02717, 0.03778, 0.04850, 0.05932, 0.07026, 0.08131, 
0.09247, and 0.10375 x 10~ 18 W m~ 2 A -1 . The filled ellipse at the upper right corner of panel / 
shows the resolutions of 0"4 in space and 60 km s _1 in velocity. The dotted horizontal line traces 
the peak velocity of the southern jet. 
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0.01 AU 0.1 AU 1 AU 10AU 100AU 

Fig. 6. — A schematic diagram showing the ranges of launching radius for a magnetocentrifugally 
accelerated outflows from a Keplerian-rotating accretion disk. The central star is assumed to have 
a mass of 1 M . The Alfven radius, the locus of which is shown in white dotted-lines, is simply 
assumed to be three times that of the launching point radius anywhere on the disk. The line of 
sight is 45° with respect to the polar axis, and the observed radial velocities (Kad) are 1/ V2 of the 
actual velocities (V). A jet with V rad =280-320 km s _1 (V=400-450 km s _1 ) is launched in a radial 
range of 0.04 to 0.05 AU, while a jet with y rad =30-200 km s" 1 (V=40-280 km s" 1 ) is launched 
in a radial range of 0.1 to 4.5 AU. The white arrow indicates the direction of the mass accretion in 
the disk. 



